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Abstract: The detrimental effect of spontaneous emission on the performance of control 
schemes designed to achieve population inversions between the ground state and a highly excited 
atomic state is studied using computer simulations. 
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We consider a A'^-level quantum system with energy levels E„ and corresponding energy eigenstates |n) for 
n = 1, 2, . . . , TV. Assuming that the system is initially in the ground state |1), our objective is to maximize 
the population difference between the ground state and the excited state \N) by applying a sequence of 
simple control pulses. For simplicity, we shall restrict our attention here to control schemes involving only 
transitions between adjacent energy levels in the model and control pulses that arc resonant with one of the 
allowed transition frequencies. It will be assumed that the allowed transition frequencies /i„ = En+i — En 
for n = 1, 2, . . . , — 1 arc distinct and sufficiently separated so that off-rcsonant effects can be neglected. 

If we neglect the finite lifetimes of the excited states and assume that each control pulse interacts only with 
its resonant transition, then it can easily be shown, e.g., using Lie group decompositions [I], that a complete 
population transfer from the ground state to the highest excited state |A^) can be achieved by applying 
a sequence of — 1 control pulses such that the nth pulse is resonant with the nth transition frequency 
Hn and has a total pulse area of 7r/(2(i„), where (i„ is the absorption oscillator strength of the transition 
n) — |n + 1). The length and shape of the pulses do not influence the outcome of the control process as long 
as the pulse areas remain fixed. For a real system, however, the decay of the excited states by spontaneous 
emission reduces the populations of the excited states and ultimately leads to the rcpopulation of the ground 
state, which is detrimental to achieving population inversion. Our aim is to study the effect of spontaneous 
emission on the outcome of the control process using computer simulations. Due to space constraints, we 
only consider a four-level model of Rubidium with energy levels and transitions as indicated in Fig. 1 as a 
concrete example. If dissipative effects are neglected then it can easily be seen that a complete population 
transfer from the ground state |I) = |5S'i/2) to the excited state |4) = I6P3/2) can be achieved by applying a 
sequence of three control pulses of frequencies /zi, /X2 and fis and pulse areas 7r/(2(ii), ■n/(2d2) and T:/{2d^), 
respectively, and that the result of the control process is independent of the shape and length of the control 
pulses. If dissipative effects such as the finite lifetimes of the excited states are taken into account, however, 
then a complete population transfer can usually not be achieved, and the performance of the control process 
depends significantly on the shape and length of the control pulses employed. 

Fig. 2 shows the difference between the populations of states |1) and |4), ^44 — pn, at the conclusion of the 
control process as a function of the total control time Tf for various ratios of the pulse lengths, as obtained 
from a computer simulation of the control process. Note that for very short control times the final yield is 
around 90% of the theoretical maximum value of one. As the total control time increases, the final yield 
P44 — pii decreases since dissipative effects accumulate. However, the computer simulation suggests that the 
final yield also depends considerably on the ratio of the pulse lengths, not just the total control time Tf. For 
our four-level Rubidium system, the naive choice of equal pulse lengths, i.e., pulse ratios I : I : f , consistently 
results in the worst yield, independent of the total control time Tf. Another obvious choice of the pulse lengths 
according to the ratio of the lifetimes of the excited states results only in a marginal improvement. However, 
the final yield and the threshold value of Tf for population inversion increase considerably as the ratio of the 
first two pulses versus the last pulse increases. These at first quite startling results can be partly explained 
by the vastly different lifetimes of the excited states of Rubidium. The lifetime of state |2} = I5P3/2) is less 
than a third of that of |3) = \^D^/2) and almost one fourth of that of |4) = I6P3/2). Hence, state |2) is the 
weakest link in the chain and it is therefore imperative to minimize the time the system spends in state 
2) in order to reduce dissipative losses. However, choosing the pulse lengths according to the ratio of the 
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Fig. 1. Energy-level diagram 
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Fig. 2. Final value of P44 — pii as a function of Tf 




10 15 20 25 

Time fnsl 




Fig. 3. Control pulse sequence (left) and corresponding evolution of the energy-level populations for an ideal 
non-dissipative system (center) and a realistic dissipative system (right) 

lifetimes of the excited states is not sufficient since the amount of time the system spends in state \n) in the 
control scheme under consideration is roughly proportional to (At„_i -I- At„)/2, not simply At„. Wc see that 
it is advantageous to choose the pulse lengths At„, n = 1,2,3, such that Ai2 + A^a > 3(Ati -f- At2), which 
explains why increasing the ratio of the pulse length of the third pulse compared to that of the first two 
pulses tends to reduce the rcpopulation of the ground state due to spontaneous emission and thus improves 
the performance of the control scheme. 

Fig. 3 shows the (envelopes of the) control pulses and the corresponding evolution of the energy-level pop- 
ulations if dissipative effects due to spontaneous emission are neglected (center) and if realistic lifetimes of 
the excited states are assumed (right). The total pulse length in the example is 30 ns, where the first and 
second pulse are both 6 ns long while the third pulse is 18 ns. The cumulative effect of spontaneous emission 
is clearly noticeable by comparing the graphs. 

The example studied in this paper suggests that it is possible to achieve substantial population inversions 
between the ground state and a highly excited state of an atom (or possibly ion or molecule) by applying a 
sequence of short control pulses, even if the finite lifetimes of the intermediate excited states are taken into 
account. In spite of inevitable population losses due to spontaneous emission, pulse lengths and ratios can 
be adjusted to minimize dissipative effects and achieve substantial yields. We believe such a scheme could be 
used to calculate the pulse configurations required to build short wavelength pulse lasers in the UV or even 
X-ray regime by applying a sequence of short optical-frequency laser pulses to a cloud of trapped atoms or 
ions to achieve population inversion between the ground state and a highly excited state. 
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